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*H . Abstract 

We employ the superpotential technique for the reconstruction of cosmological models with 
a non-minimally coupled scalar field evolving on a spatially flat Friedmann-Robertson- Walker 
■ background. The key point in this method is that the Hubble parameter is considered as a function 

CN| , of the scalar field and this allows one to reconstruct the scalar field potential and determine the 

dynamics of the field itself, without a priori fixing the Hubble parameter as a function of time 
or of the scale factor. The scalar field potentials which lead to de Sitter or asymptotic de Sitter 
solutions, and those which reproduce the cosmological evolution given by Einstein-Hilbert action 
plus a barotropic perfect fluid, have been obtained. 
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1 Introduction 



There are two possible ways of describing the accelerating expansion of the Universe [I]. One can 
assume that General Relativity is the correct theory of gravity and introduce a smoothly distributed, 
slowly varying cosmic fluid with negative pressure called dark energy [2|. Alternatively, one can 
modify the theory of gravity [31 El [5]. Some of these modified gravity models, for example, F(R) 
gravity models, can then be mapped into general relativity with additional scalar fields by a suitable 
conformal transformation of the metric (see e.g. [5j [6]). It is well-known that scalar fields play an 
essential role in modern cosmology since they are possible candidates for the role of the inflaton 
field driving inflation in the early universe [8] and of the dark energy substance [2] . The universe 
expansion history, unifying early-time inflation and late-time acceleration, can be realized in scalar- 
tensor gravity [9]. 

The modified gravity models, as well as cosmological models with scalar fields, often include some 
function which cannot be deduced from the fundamental theory. It is then a natural question to ask, 
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why some specific form of such a function is chosen and what is the physical motivation for this choice. 
The technique of reconstruction has attracted the attention of researchers for a long time and has been 
developed for different cosmological models QH]-[3S]. For example, the reconstruction of potentials 
for models with minimally coupled scalar fields has been studied in |10]-[20j and, in the two-field 
case, in |21[ |2"2"] 123] , while a similar procedure for tachyon models was discussed in |24[ |2"5"| l2Uj 127] . 
The reconstruction of potentials for scalar fields non-minimally coupled to gravity was considered in 
[28, 29, 30j. We can also mention the reconstruction procedures in models with non-minimally coupled 
Yang-Mills fields [31], in F(R) and Gauss-Bonnet gravity models \20\ I32j. in F(T) models, with T 
being the torsion scalar [33], in the nonlocal gravity model |34j and its local formulation |35| . 

Gravity models with non-minimally coupled scalar fields, are actively studied in cosmology |36j- 
[43] (see also [5] and references therein). In particular the models with the Hilbert-Einstein term plus 
the term proportional to the Ricci scalar multiplied by the square of the scalar field were intensively 
studied in inflationary [33] and in quantum [45] cosmology. 

Generally these models are described by the following action 



where U(a) and V(a) are differentiable functions of the scalar field a. We use the signature (+, —,—,—) 
and g is the determinant of the metric tensor g^ u . 

The reconstruction procedure for the induced gravity models (U(a) = £cr 2 where £ is a constant) 
has been proposed in [29]. In such a case it has been shown that one can linearize all the differential 
equations which should be solved in the reconstruction procedure to get the potential corresponding to 
a given cosmological evolution. This property allows one to obtain the explicit forms of potentials re- 
producing the dynamics of a flat Friedmann-Robertson-Walker (FRW) universe, driven by barotropic 
perfect fluids, by a Chaplygin gas and by a modified Chaplygin gas [2"9"] . 

In this paper, we consider another reconstruction procedure for the models described by the above 
action ([1]). Such a method is similar to the Hamilton-Jacobi method (also known as the superpotential 
method) and is commonly applied to cosmological models with scalar fields \19\ [2T| [22], [23] [27] as well 
as in brane models [56]. The key point in this method is that the Hubble parameter is considered as 
a function of the scalar field a. 

The two methods described above supplement each other and together allow one to construct dif- 
ferent cosmological models with some required properties. In particular, the reconstruction procedure 
illustrated in [29] is useful in order to obtain the potential and the explicit evolution of the scalar 
field when the Hubble parameter is given explicitly as a function of the cosmic time t or of the scale 
factor a. The superpotential method also allows one to do so (see Section 4), but for U(a) = £cr 2 
the procedure proposed in [29] is simpler. At the same time, however, the superpotential variant of 
the reconstruction procedure allows one to construct models and their exact solutions without any 
knowledge of the exact behavior of the Hubble parameter. For example, in Section 6 we describe the 
dynamics of a model with the Hubble parameter evolving to a nonzero constant value at late times. 
Such a method also allows one to calculate the function U(a), when the Hubble parameter is given 
as function of a and the behavior of a as a function of the cosmic time t (or the scale factor a) is 
known. Specifically, we assume that 4? = F(<j), where the function F is given. So, the superpoten- 
tial procedure allows one to reconstruct not only the potential, but also either U(a), or the Hubble 
parameter, if the evolution of a is given. Furthermore we can use this method to reconstruct scalar 
field potentials with some given property. For example, in Section 7 we find the conditions on the 
functions U (a) and F(a) for which the potential obtained is polynomial. 

By using this method we find new exact solutions in the induced gravity case. We also applied 
the procedure to cosmological models with a non-minimally coupled scalar field described by ([I]) and 
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U(a) = £<t 2 + J, where J is a constant. In this physically important case, the method, proposed in [29] . 
requires one to find solutions of second order nonlinear differential equations, whereas on using the 
superpotential method one only has to solve linear differential equations and first order autonomous 
nonlinear differential equations. Thus, the solutions can be always obtained at least in quadratures. 
One important case treated in detail is the reconstruction of de Sitter solutions. In such a case, the 
method proposed in [29] is still useful for arbitrary U(a) (see Section 3). 

The paper is organized as follows. In Section 2, we review the basic equations for a gravity model 
with a non-minimally coupled scalar field in a FRW background and we describe the algorithm of 
the superpotential reconstruction procedure. In Section 3, we study the de Sitter solutions by using 
the methods proposed in [29] and in this paper and compare results. In Section 4, we find the non- 
minimally coupled models reproducing the dynamics of a flat FRW universe, driven by barotropic 
perfect fluids. In Section 5, we assume that the scalar field is proportional to the hyperbolic tangent 
(or cotangent) of the cosmic time and consider possible behaviors of the Hubble parameter and the 
corresponding potentials. In Section 6, we consider the induced gravity model and find solutions and 
potentials with a Hubble parameter which evolves towards a constant value at late times. In Section 7, 
we find conditions on the functions U(a) and F{a) sufficient to obtain a model with a polynomial 
potential. Finally Section 8 is devoted to the conclusions. 

2 Cosmological models with non-minimally coupled scalar field 
2.1 Homogeneous Equations 

Consider the evolution of an homogeneous scalar field on a spatially flat FRW universe with the metric 

ds 2 = dt 2 - a 2 (t) (dxj + dx\ + dxf) . 
The Einstein equations derived from the variation of the action ([I]) have the following forrrfl |29j : 

6UH 2 + 6UH = hr 2 + V, (2) 

2U (2H + 3# 2 ) + AUH + 2U + ^& 2 - V = 0, (3) 

where differentiation with respect to time t is denoted by a dot, and the Hubble parameter is the 
logarithmic derivative of the scale factor: H = a/a. The variation of action ([I]) with respect to a gives 
the Klein-Gordon equation: 

a + 3H& + V a = &[H + 2H 2 ^j U >(7 , (4) 

where the subscript >cr indicates the derivative w.r.t. the scalar field a. Combining Eqs. ([2]) and ([3]) 
we obtain: 

iUH - 2UH + 2U + a 2 = 0. (5) 

This equation plays a key role in the reconstruction procedure. 

1 In [29], the authors have used a convenient choice of Newton's constant G, 8nG/3 = 1, and dimensionless quantities. 
Alternatively, one can assume that G is a part of U(a) and do not fix the value of this constant. 
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2.2 The Hubble parameter as a function (superpotential) of the scalar field. 

Let H = Y(a), then (0) takes the following form: 

iUY a a - 2YU a a + 2[/ CT(T cr 2 + 2U a a + a 2 = 0. (6) 

Let us introduce the function F(a) defined as 

a = F(a). (7) 

Substituting a and a = F^F into ©, we obtain: 

4UY a + 2(F t<T - Y)U a + (2U aa + 1)F = 0. (8) 

Equation (JSj) contains three functions. If two of them are given, then the third one can be found as the 
solution of a linear differential equation. Let us note that Eq. (J8|) is a first order differential equation 
both for Y and F . 

If, for example, U and F are given, then Eq. ([8]) can be easily integrated to obtain Y{a): 



Y(a) 



d ° ^372 + Co ) VU 



(9) 



where cq is an integration constant. In a similar fashion, for given Y(o) and U(cr), we obtain F{a) by 
integrating (jSJ): 



F(a) 



da 



Y{a) - 2U(a) 



dU(&) 
da 



+ c e" 



-T(a) 



(10) 



where 



20 + 1 
2 1 

da 



da 



and Co is an integration constant. 

The potential V{a) can then be obtained from 



V{a) = 6UY 2 + 6U, a FY - -F 2 . 



(11) 



In order to find the time evolution a(t) we finally integrate Eq. ((7|), which can always be solved in 
quadratures. 

The evolution of the Hubble parameter H and of the scalar field a can be expressed either through 
the cosmic time t, or the scale factor a. In the latter case it is more convenient to adopt the dimen- 
sionless parameter N = ln(a/ao) as independent variable. In terms of N one has to re-express the 
time derivatives as 



d 2 

and \=H 2 



r/- 



1 



./ d 



dt " dN — " dt 2 " dN 2 2 1 ; diV 



(12) 



where the prime denotes the derivative with respect to N. Equation ([8]) then takes the following form 



2U{H 2 )' - 2U'H 2 + 2U"H 2 + (H 2 )'U' + {a'fH 



II tt2 



t2\IttI 



J\2 tt2 



0. 



(13) 
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On denning the functions W(o) and F(a) as follows 

W(a) = H 2 

and 



we find 



F{a) = a', 



da dN 

U' = U a F, U" = (u <aa F + U^F a ) F 
and Eq. (|13p can then be cast in the form: 



2UW, a - 2U a W + 2 (U, aa F + U a F, a \ W + W^U^F + FW = 0. 
If, for example, F(a) and U(a) are given, the function W(a) can be determined by integrating: 



(14) 

(15) 

(16) 
(17) 

(18) 



w, a _ 


2U a - 


2 (u >aa F + - F 


w 




hu + U a F\ 



(19) 



and finally the potential can be reconstructed by using 

V = W 



6U + 6U a F - -F 2 



(20) 



Equation (fl~9j) connects the functions W, U, and F and is equivalent to Eq. ([8]). Indeed, on substituting 
F = YF into Eq. © we get 



2^ 



2U + U a F 



(21) 



The latter equation is equivalent to Eq. (fl~9|h because 1^(0") = ^ 2 (o"). We thus observe that the 
superpotential method allows one to unify the reconstruction procedures for the functions in ([8]) given 
as functions of t or of N. Indeed, one can start from Eq. (|8|) and find suitable functions Y, F and U 
satisfying it and, after that, it is trivial to get W and F. 



2.3 Quadratic non-minimally coupled models 

In this paper we consider 

U(a) =t<j 2 + C 1 a + J . 



(22) 



The case of induced gravity, for which the reconstruction procedure has been proposed in [29] , corre- 
sponds to a particular choice of U(a) given by C\ = and Jo = 0. 
Note that the quadratic polynomial (I22p can be rewritten as 



U(a) = £(o" — ctq) 2 + J, with ao 



J = Jo 



4£ 



(23) 
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In this paper, we consider the case <to = (equivalent to C\ = 0). Let us observe, however, that such 
a restriction does not lead to any loss of generality because any solution found with oq = can be 
extended to the <7o / case simply through the change of variable a — > a — oq. 
For U(a), given by ([22]) with uq = 0, we obtain from Q and (fTUj) : 



Y(a) 



A£aF d + (4£ + 1)F 



4(£a 2 + J)3/2 



Ve^ 2 + J, 



F(a) 



( <r 2 + | ) 



a 



-1-1/(4?) 



(24) 



(25) 



Let us note that the behavior of the Hubble parameter can be very different in a model with zero 
and nonzero J. To illustrate this point let us consider the example with an arbitrary linear function 
F(a): 

F = fia + /o, (26) 
where are constants. From Eq. (|24h . we get: 



(27) 
(28) 



Yj(a) = £ + 2h - (4C ±^ a + cqv^TJ , J / 0, 

v / \ /l , of . ( 4 £ + l)/o . 7 n 

F °(°") = 77 + + 77 + c cr , J = 0. 

Therefore, if /o 7^ and £ 7^ —1/4, then Yo(°") is n °t a limit of Yj(a) at J — >• 0. 



3 Models with de Sitter solutions 

In this section we shall consider the general form of the potential V(<r), which leads to the existence 
of the de Sitter solution 

H = H = const (29) 

for an arbitrary non-minimal coupling U(a). In this case, it is convenient to apply the technique 
developed in [2T?] . 

First of all let us notice, that if the scalar field a does not depend on time t (or, in other words, 
on the scale factor a) then all the considerations presented in the end of the second section of paper 
[29J are still valid and we obtain immediately the de Sitter solution for the model with the potential 
(see Eq. (31) of the above cited paper) 

V{a) = V U 2 (a), (30) 

where Vo is an arbitrary positive constant. The constant value of the Hubble parameter is then (see 
Eq. ©) 

For the case of a time-dependent scalar field a, we use Eq. (19) from paper [29] and we rewrite it 
in terms of N = lna/ao 

o-" J 1 ~' H ' 1 t„i\2 ( l + 2 lb^ \ 1 2UH ' n { on\ 

V da J da 
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For the de Sitter case, this equation reduces to 



a" - a 1 



1 + 2 



d?U 
da' 2 



ydU 
' da 



{a'f- 



(33) 



Using ^ = a after some algebra, one can integrate both sides of Eq. (|33j) and obtain: 



lno-'-TV 



da . fdU\ 

da v 



(34) 



where -Bo is an integration constant. 

We now exponentiate and integrate the above equation to find exp(iV) as a function of the scalar 
field a: 



exp(iV) = j-Jda^exp [\ J d *w\ + 



(35) 



where Bq = exp(Bo) and B\ is an arbitrary constant. 
Consequently we obtain a' as a function of a: 



a 



ex P ( -iJdvik 

da 



dU 
da 



BiB + I da*^exp | i J da^j 



(36) 



and the corresponding potential can be found by substituting (|36j) and (|29j) into Eq. (1201) . 
For the particular case: 

U(a) = £a 2 + J, 



we obtain 



Substituting Q38} into © we get 



4£ Bi^o 

:<r H ~ — a 4 « 



1 + 8£ 2£ 



(37) 
(38) 



v = m 



gJ 2(3 + 32Q(1 + 12Qg j2 4^0(1 + 120 ^.^ _ B(B^_*J+1 



(8C + 1) 2 



l + 8£ 



(39) 



Let us apply the superpotential method to the same model with the de Sitter evolution. Using ([10 
we obtain 



F(a) 



and, in particular, for U, given by (|37|) . 



Hoe T da + cq 



(40) 



(41) 



Therefore, 



a 1 = F(a) 



F(a) = 4£ jp_ -i-i/(40 
y(<r) 8£ + l #o 



(42) 
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This formula coincides with (|38p . if we choose Co = HqB\Bq/(2£), so we conclude that the potential 
can be reconstructed with both the procedures. 
Solving Eq. ([7D, we obtain 



a(t) 



a Q e^ + £2! 



8 ^ + l)l^/(8€+D 



Hot 



(43) 



where do is an arbitrary constant. 

For the case cq = 0, F(p) is a linear function. The potential V is the following quadratic polyno- 
mial: 

t2 



v = 2m 



3J | (3 + 320(1 + 12^ ^ 



(se + i) 2 

For £ = —1/12 and £ = —3/32 the potential V is a constant. 

The case £ = —1/8 should be considered separately. In this case, using Eq. (fTUj) . we get: 

F(a) = aH ln ( — 



(44) 



where ctq is an integration constant, and the corresponding potential has the following form: 



(45) 



24 J - — I In ( — ) +3 + v / 3][ln( — ) +3-^3 



o-o/ 



The scalar field evolution is given by 



a(t) = do exp 



D H (t-t ) 



(46) 



It is easy to see that in the case of induced gravity (J = 0) the expressions (f3Uj) and (|4"4"|) can be 
derived from the expression (51) in the paper provided we identify c\ = — BlB °^ +8 & anc j ^ _ 2£. 
Also the particular cases 7 = —1/6 (conformal coupling) and 7 = —1/4 corresponding to £ = —1/12 
and £ = —1/8 respectively were considered in |29j . 



4 Power-law Solutions 

In general relativity with a barotropic perfect fluid, the Hubble parameter evolves as 



-3(tu+l) 

fl a (a) = flg(-l =^ 2 e nJV , (47) 



a 

where the constant w is the equation-of-state parameter of the fluid and n = —3(w + l). Let us note 
that the Hubble parameter as a function of time is H(t) = — ^ = 3 ( w + 1 ) f ■ Power-law solutions for 
induced gravity models (J = 0) and the corresponding potentials have been studied in [29]. Let us 
note that, in the case with J ^ 0, the method employed in [29] leads to a differential equation which 
cannot be solved analytically. For such a case it is convenient to adopt the superpotential technique 
for reconstruction. 

Let us start from Eq. (|19p and observe that the power-law dependence of H 2 on the scale factor 
gives n = (H 2 )' / H 2 , or equivalently 

-=¥■ <«> 
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Then Eq. (|T9|) becomes 



n 



2U a - 2 U„ a F + U~F a 



(49) 



F [2U + U a F } 

After substituting (|37p . Eq. (|49p can be recast in the following form: 

4£ctF ct F + 2£a{n - 2)F + (4f + l)F 2 + 2nJ + 2n^a 2 = 0. (50) 

In particular, for n = 2, Eq. (|5Q|) is a linear differential equation for F 2 which can be solved exactly 
leading to the following solution: 



For this solution from Eq. (148 P. we get 



W((r) = VF e J . (52) 



For the values of £ leading to a singular expression (|5ip . Eq. (|5Q|) has the following general solution: 
• for £ = -1/4 



F(a) = ±j8J\n( — )-a 2 , (53) 



for £ = -1/8 



F(o-) = ±J-2a 2 ln(— ) - 8J (54) 



where do is the integration constant. 

It is impossible to solve Eq. (|50p by the separation of a and F for the case of an arbitrary n. Hence 
we look for particular solutions of Eq. (|50p by making a suitable ansatz for F as a function of a and 
then algebraically fixing the free coefficients of the model. We choose the ansatz 

F(a) = A + A x a , (55) 

by which Eq. (150p reduces to a quadratic polynomial of a, which can be solved by requiring that the 
coefficients of the diverse powers of a are zero. Let us note that the above ansatz (|55p when A\ ^ 
corresponds to the following evolution for the scalar field in terms of N 

a(N) = B exp (A±N) - -± (56) 

A\ 

where B is an integration constant. It is then straightforward to obtain the evolution a(t) by inte- 
grating: 

dN 2 f t \ ~ 2/n 

dN = —dt = Hdt = dt N = In — . (57) 



dt nt \to, 

and substituting into Eq. (|56p . Equation (|50p gives three constraints on the 5 parameters of the 
model: Aq, A±, n, £ and J. If we assume that the function U(a) is given by (f3Tj) . or in other words 



9 



that the parameters £ and J are fixed, we then find the following expressions for the remaining three 
parameters: 



-4 



A, 



-14f - 44£ 2 - 1 ± (6£ + 1)V 52 ^ 2 + 16£ + 1 



J 



(4£ + l)2£ 



-36g 2 - 12g - 1 ± (6g + l)\/52£ 2 + 16£ + 1 
(6£ + l)(4£ + l) : 



(58) 



n 



Mg + 44£ 2 + 1 T (6g + W52g + 16^ + 1 
(4£ + 1)£ 



These solutions exist only if £ ^ —1/4 and £ ^ —1/6. The parameters n and A\ are real if and 
only if £ < (-4 - V3)/26 or £ > (-4 + y/3)/26. Let us note that A = if and only if J = 0. 

On the other hand, if we assume, that n is given and seek both the functions F and U then Eq. 
(|50p leads to two nontrivial solutions with unconstrained and: 



A 2 - 



6,2 



-An - a ± (n - 2) yfa 
An (a — 2n) 

8n — a =p (n — 2) -^/a 



Ai 



1,2 



8 (a - 2n) 
2 — n =F 



(59) 



(60) 



(61) 



with a = A — 20n + n 2 and where Aq and — ^4o correspond to the same values for the other parameters. 
These solutions describe a Universe evolving according to (|47p due to the presence of a non-minimally 
coupled scalar field with the coupling to the Ricci scalar given by (1601) . Note that a is non negative 
and expressions (f59"|) - (f6"Tj) are real provided n > 10 + A\/6 ~ 19.8 or n < 10 - 4^/6 ~ 0.2. Negative 
values of n correspond to a decreasing H in an expanding universe and describe relevant cosmological 
evolutions. The limit n — > gives the de Sitter solution and the solutions with small, positive n still 
describe interesting cosmological evolutions corresponding to a super accelerated era (the fact that the 
Universe is currently undergoing a phase of super accelerated expansion is not excluded by supernovae 
observations |48j). 

Equation (j50|) also has 4 more solutions corresponding either to the induced gravity or the general 
relativity framework. We find induced gravity solutions (with J = 0) having Aq = 0, £ unconstrained 
and 



At 



3,4 



-(n-2)e± v / e(«e-2n) 



(62) 



l + 8£ 

or £5 = —1/8 and Ai-% = n/ (2 — n). In these cases the scalar field evolution is given by cr(N) = 
B exp A\N. The general relativity solution is finally for A\ = 0, £ = 0, J = — j4g/(2n) and a scalar 
field evolution given by a = AqN + B where B is an arbitrary constant. The constant Aq is physically 
constrained by the value of the Planck mass in the Einstein-Hilbert action. 
Eq. ([35)1 can be solved for W with F given by the ansatz ([55]) : 



W = Wo exp 
where Wo is an integration constant 









[/>: 


= Wq exp 


[J 







n 



A + Aid- 



da 



W (A + A 1( j) a i . 



(63) 
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On using Eq. (|20|) . we obtain the corresponding potential 



V = W {A + Aia) A i 



6 (J + £<r 2 ) + 12£cr (Ao + A x o-) 



(A) + A x a) 



(64) 



with J, £ and given by (f59l) - (f6T|) respectively. In the induced gravity case (solutions 3 and 4) the 
potential takes the following form: 



n , 2 _ (6-n)e± v /g( Q ;-2ra) 

y 34 = VFoa Al ;3. 4 =Wb<r 2 « (65) 

where we absorbed all the parameters and the integration constant Wq into Wq. Correspondingly the 
evolution of the scalar field is given by 



03,4(^0 = -Bexp 



-(n-2)g± v^(qg-2n) 

i + 8e 



N 



(66) 



These two solution where already found in [29] . In the n — > (and a — > 4) limit one then obtains the 
de Sitter solutions (see [2S]) with 



(dS) 



W a z 



and 



4 dS \N) = Bexp 



-A/ 



= W O 

a f 5) (iV) = B. 



(67) 
(68) 



1 + 8^ 

If we consider the induced gravity solution left, with £5 = —1/8, we are led to the potential 

V 5 = W^- n (69) 
and the evolution of the scalar field is given by 

n 

Let us note that such a solution cannot be obtained from the solutions 1-4 with an appropriate limit 
procedure. 

The general relativity solution has 



a 5 (N) = Bexp 



-N 



(70) 



W = Wq exp 
and a corresponding potential of the form 

V = Wq exp 



n 

A 



-a 



a 



A 



-a 



(71) 



(72) 



which is the well-know exponential potential associated with power-law expansion. 

The solutions 1 and 2 deserve some more discussion. In the limit J± ; 2 oc Aq — > one obtains 
induced gravity with the non- minimal coupling fixed at £ = £1^. In such a case the potential takes 
the form 

V 1>2 = W a 4 (73) 

and the scalar fields evolves as 

2 — n =F y/a 



(?i,2(N) = B exp 



(74) 



Note that the solution 1, in the induced gravity limit, is the solution 4 with £ = £1 and the solution 2, 
in the same limit, is the solution 3 with £ = £2- Thus, when J = 0, the solutions 1 and 2 are particular 
cases of the solutions 3 and 4. 
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5 Solutions with the hyperbolic tangent 

Let us construct cosmological models, when the Hubble parameter is a function of the hyperbolic 
tangent. Such solutions are popular in cosmology, because they can describe a bounce [U] and late- 
time acceleration [171 [21] having de Sitter /anti-de Sitter attr actors both in the past and in the future. 
Such models can be reconstructed, for example, starting with the following ansatz for the scalar field 
evolution 

a(t) = A tanh [u(t -t )], (75) 

where A, oj and to are constants. Note that to can be complex, so the parametrization (175j) includes 
the functions a(t) = ^4coth [co(t — to)] as well. For such functions Eq. (|7|) takes the form 

F(a) = u (a- -^ 2 ) . (76) 

If U(a) = £cr 2 , then after substituting F(a), given by (f76|) . into formula (f24"|) we can integrate the 
resulting equation to obtain: 

where do is an arbitrary constant. The function Y(a) has the simplest form for £ = —1/12: 

Y(a) = c a-—. (78) 
a 



^JsqAjOV. (79) 



The corresponding potential is the fourth degree monomial: 

2A 2 ~ 

For Co = 0, we obtain the following evolution for the Hubble parameter: 

H= -wcoth[w(t-to)] (80) 

and conversely we get H = — wtanh[u;(t — to)], when a(t) = ^4coth[w(t — to)]. 
In the more general case with U(a) = —a 2 /12 + J we get: 

y = c oV / 12J-a 2 -^. (81) 
6J 

For Co = 0, the Hubble parameter is proportional to a and the potential has following polynomial 
structure: 

(6J + A 2 ) 2 C0 2 4 [A 2 \ 2 2 1 2 ,„ 

Another way to get the desired Hubble parameter evolution (with constant attractors in the 
past and in the future), is to express H as a function of a evolving as (|75p and then solve the Eq. 
as a differential equation for U{a). On assuming 

iJ = Y(a) = B-Ca, (83) 

where -B and C are constants, and then substituting (fTBT) and (i83l) into Eq. (jHJ), we get the following 
second order linear differential equation for U(a): 

2n(A 2 - a 2 )U :aa + 2 [(C - 20)cj - S] f7, CT - 4CC7 + (^ 2 - <j 2 )tt = 0, (84) 
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where we set u = OA 

A particular solution of this equation is the second degree polynomial 



U{a) 



1 



2 , 

12 6(20 + C) 



B 2A 2 CQ + 4A 2 n 2 - B 2 

CT+- 



12(20 + C)C 



(85) 



having the form (|22H . For any values of A, B, C and 0, except = — C, the potential V(a) is a fourth 
degree polynomial with a coefficient of the high-order degree monomial equal to — (0 + C) 2 /2. 
Let us consider, for example, the particular set of parameters B = A and (7=1. We get 



1 



U{a) =--2 a2 + 6(20 + 1) 



A (20 + 40 2 - lU 2 



12(20 + 1) 



(86) 



hence, U(a) = at 



0-12 = ^4 



1 ± 2V20 2 + 20 3 



20 + 1 

It is easy to see that U(A) > for all O > 0, so, if we choose as a solution a(t) = Atanh(wt), then 
U(t) is positive at late times if A > 0. We further observe that (see Fig. [T]), when = 1, U(a(t)) > 
at any time because —A ^ a(t) ^ A. For O = 1, we get 




-1,5 



-0,15 




Figure 1: The function U(a), given by (i86"1) . (left) and the potential V(a), given by (|87|) . (right) at 
A = 1, B = 1, (7 = 1 and O = 1. 



VYcr) = - 2a 4 + -Aa 3 + -^a 2 - -A 3 a + -^l 4 
V ' 3 3 3 3 



(a + A)(3a- A){a - Af 



(87) 



and, by the change of variable a = a — A/3, the expression (|86|) becomes 

U(cr) = - -a 2 + -A 2 , 
{ J 12 27 

and takes the form of (|3T[) with ^ = — 1/12. In terms of a, we finally get 



(89) 
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So, we found a model with exact solutions and U(a) in the form U(a) = £<r 2 + J. We note that 
if we consider general values of A, C, and u still the particular solution (|85p is for £ = —1/12. This 
value of £ corresponds to the case of the conformal coupling |29j. 



6 Non-monotonic behavior of the Hubble Parameter in the Case of 
Induced Gravity 

In the previous sections we considered the case of U(a), specified by (|57|) . The case of induced gravity 
(J = 0) has been considered in detail in [29], where another method of reconstruction has been used. 
In this section, we find new induced gravity cosmological models with exact solutions on using the 
superpotential method. We demonstrate that the same evolution a(t) leads to exactly solvable models 
with different potentials and different qualitative behavior of the Hubble parameter. 
Given U{&) = £cr 2 , and taken Y(a) as a generic polynomial: 



v 



Y(a) = ^A k a k , 



(90) 



k=0 



where A k are constants, from (|25j) one can obtain F(a) and show that it does not depend on A\. For 
N = 2, we obtain 



_ 4g (flgg + 1)4> - gg + l)A 2 a*) a 
{) (8£ + l)(16£ + l) +C ° a 



(91) 



When Co = 0, F(a) is a cubic polynomial and Eq. ([7|) has the general solution which can be written 
in terms of elementary functions: 



a{t) 



yXjg + l)Mgc 2 e-" f + (8£ + l)(16g + gggc 
(16f + l)A)C 2 e-^ + (8£ + 1)A 2 



(92) 



where w = 8£^4o/(8£ + 1), c 2 is an arbitrary integration constant and we assume that £ ^ —1/8. 

The corresponding potential, V(cr), is the sixth degree polynomial which, for example, when £ = 1 
has the following form: 



__, N 910 156 ,, 5 

V(a) = Aia b H AiA 2 a 

y ' 289 2 17 



QA\ 



2236 
153 



A A 2 



^4 , 52 3 910 2 , 



Let us analyze the cosmological consequences of the solution obtained. If oj > 0, then 



lim a(t) 

t— voo 



(16g + l)y^ 
(8£ + l)V^ 



(93) 



In the case oj < 0, the function a(t) tends to zero at late times. Hence, the Hubble parameter tends 
to a constant at late times for any casdE On Fig. [2] we consider H(t) for different values of A\, for a 
positive, £ = 1, A 2 = 1, and ^4o = 2. The solution a(t) in ([92]) is associated with different behaviors 
of the Hubble parameter. In particular, at A\ = —4 and A\ = —2 we get a non-monotonic behavior 
of H(t). The behavior of the Hubble parameter, similar to the case A\ = —4, has been found in 
the quintom models with the sixth degree polynomial potential |21j . Such behavior can, in principle, 



2 We choose Y(a) to be a quadratic polynomial and therefore A2 7^ 0. 
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Figure 2: The function H(t), given by ([50]). with Ai = -6, A\ = -4, A\ = -2, and Ai = (from left 
to right). At all pictures we use A2 = 1, Aq = 2, and C2 = 100000. 

describe both inflation and the late-times accelerating expansion of the Universe and is thus physically 
relevant. 

In this section, we have obtained new induced gravity models, with a polynomial potential and 
physically interesting behavior of the Hubble parameter. Let us note that the explicit form of the 
Hubble parameter H(t) is too complicated to be guessed and used in the reconstruction approach, 
proposed in [29]. We have obtained a three-parameter set of the induced gravity models with sixth 
degree polynomial potentials. The freedom of choice of parameters A{ allows one to set additional 
conditions on the model. 

7 Construction of models with polynomial potentials 

On considering cosmological models with scalar fields, which are inspired by some fundamental theory, 
for example, string field theory, it is difficult to get the exact form of the potential, but it is possible 
to get, at least, some properties of the potential. In (TTl [21] string field theory inspired models with 
minimally coupled scalar fields and polynomial potentials have been considered. Let us address the 
question regarding what functions F and U are associated with a polynomial potential. 

To get a polynomial potential from formula (llip we assume that F{a) is an arbitrary N degree 
polynomial: 

N 

F(a) = Y J B n v n - (94) 

fc=0 

If U = £cr 2 , then from Eq. (jSJ) we get 

Y(a) = ^ (l + ^) ~ |f(12* + ^ ln ^ + Y ^ ( 95 ) 

where Yp(a) is a polynomial. The term proportional to cr -1 gives the constant term in the potential 
(jlip . thus only the term proportional to a ln(«r) should be eliminated to obtain a polynomial potential. 
Hence, we come to the conclusion that for an arbitrary polynomial F{a) we get a polynomial potential 
if £ = —1/12. This value of £ corresponds to the case of conformal coupling |29j. For other values of 
£ we get a polynomial potential only if B2 = 0. 

If U is an arbitrary quadratic polynomial, we get the following result: the coefficients of Y(cr), 
proportional to Bq and B\ are always polynomial, those proportional to B2 and B3 are polynomial 
if £ = —1/12. The coefficient proportional to B4 includes logarithmic terms. The function Y(a) also 
contains the term, proportional to ^/U(a), multiplied by an integration constant, therefore, in order to 
get a polynomial V(a) this integration constant should be zero. Hence, to get a polynomial potential, 
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we should either set N = 1 or N = 3 and £ = — 1/12. In the latter case, we obtain a sixth degree 
potential for nonzero B%. For B3 = (or equivalently N = 2), the degree or the potential can not be 
more than four. 

8 Conclusion 

In this article we employed the reconstruction procedure for cosmological models with non-minimally 
coupled scalar fields evolving on a flat Friedmann-Robertson- Walker background. Besides their cosmo- 
logical implications, models with non-minimally coupled scalar fields, including the Hilbert-Einstein 
term, are interesting because of their connection with particle physics. In particular we refer to models 
of inflation, where the role of the inflaton is played by the Higgs field non-minimally coupled to gravity. 
Such models have great relevance and were intensively discussed during the last years [19] . 

It is difficult to apply the reconstruction procedure proposed in [29] in the induced gravity context 
to more general modified gravity models. For these models, with the superpotential method, one 
can find many exact solutions for the homogeneous scalar field-gravity system and the corresponding 
potential of the scalar field. Such a method translates the problem of integrating a second order 
differential equation with time dependent coefficients into the integration of a first order differential 
equation once either the Hubble parameter is expressed as a function of the scalar field or the behavior 
of the scalar field itself is given. Let us note that both the variants of the superpotential technique we 
are referring to do not need the expression of the Hubble parameter in terms of the cosmic time or 
of the scale factor. This method has been actively employed in cosmological models with minimally 
coupled scalar fields [191 [ZD E21 [23], but, to the best of our knowledge, it has not been applied to 
cosmological models with non-minimally scalar fields. Hence the superpotential method is a powerful 
tool to find exact analytical solutions of cosmological interest. 

In this article we used this method to find exact analytical solutions leading to some physically 
interesting behaviors of the Hubble parameter. First we have found the potentials and the correspond- 
ing evolutions of the associated scalar field leading to de Sitter and power-law solutions. Power-law 
solutions reproduce the expansion of the Universe driven by a barotropic perfect fluid. For the for- 
mer case we could formally reconstruct the dynamics of the scalar field-gravity system on assuming 
a generic coupling between the field and the Ricci scalar. In such a case both the superpotential 
technique and the reconstruction procedure proposed in [29] can be applied. We then specialized the 
results to the non-minimally coupled case with a general coupling constant £. The case of power-law 
expansion is more involved and we could only find and analyze some solutions with the parameters 
suitably tuned. Still, however, the reconstruction could only be performed easily by means of the 
superpotential technique. 

We further investigated a few models having a de Sitter/anti de Sitter attractor far in the past or 
in the future. In the non-minimal coupling case the dynamics of a scalar field with a suitable quartic 
potential has been shown to have these attractors. In the induced gravity limit, we then also found a 
set of sixth degree polynomial potentials leading to a cosmic evolution with a constant H{t) behavior 
both in the past and in the future. 

Polynomial potential frequently appear as a result of the reconstruction procedure because poly- 
nomials are integrable and they can be systematically used for obtaining exact solutions with the 
superpotential technique. Hence we finally analyzed the sufficient conditions on the functions U and 
F to get a model with a polynomial potential. 
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